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ABSTRACT 

Stars of all evolutionary phases have been found to have excess infrared emission due to the 
presence of circumstellar material. To identify such stars, we have positionally correlated the 
infrared MSX point source catalogue and the Tycho 2 optical catalogue. A near/mid infrared 
colour criteria has been developed to select infrared excess stars. The search yielded 1938 
excess stars, over half (979) have never previously been detected by IRAS. The excess stars 
were found to be young objects such as Herbig Ae/Be and Be stars, and evolved objects such 
as OH/IR and carbon stars. A number of B type excess stars were also discovered whose 
infrared colours could not be readily explained by known catalogued objects. 

Key words: - stars: circumstellar matter - early type - evolution - pre-main sequence 



1 INTRODUCTION 



The discovery of excess far-infrared emissio n from a Lyrae by the 
Infrared Astronomical Satellite (IRAS) by lAumann et alj il984) 
demonstrated for the first time that proto-planetary material around 
other stars could be both detected and studied. Further studies of 
IRAS point sources revealed that many types of star displayed 
strong emission at IRAS wavelengths. This excess is mostly due 
to thermally re-radiating circumstellar dust or, mostly in the case 
of hot Be stars, due to free-free and bound-free emission from their 
ionised gaseous disks. Con sequently, the IRAS Point Source Cat- 
alogue tBeichman et aljl988i) became a much used source for the 
search and identification of stars surrounded by circumstellar ma- 
terial. The use of an IRAS two-colour diagram to systematically 
identify evo lved stars in the IRAS PSC, was d eveloped and un- 
dertak en bv Ivan der Veen & Habina ( 198^) and IWalker & Cotier] 
il988l) . They found that the majority of the infrared (IR) excess 
stars were carbon or oxygen rich asymptoti c giant branch (A GB) 
stars undergoing mass loss. Other studies bv lPlets etaf]<1997t) and 
[jura l,19_99) also discovered giant and first ascent red giant stars 
with IR excess. 

Another manner of finding objects was to use optical cata- 
logues and cross-correl ate those with the IRAS PSC. Systematic 
studie s suc h as those by Pottasch et al. l'T98?),'stencelTBactaa3 
il99ll) and lOudmaiier et al. (.199Z) proved very successful in re- 
turning all well-known Vega-type systems and in addition uncov- 
ering large numbers of both young and evolved stars alike. For ex- 
ample, the presence of IR excess became one of the defining char- 
acteristics of pre-main sequence Herbig Ae/Be stars, and evolved 
post-AGB stars where the AGB mass loss phase had ended and 
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the IR excess t races the cool, deta ched dust shell (see for example 
the reviews bv lvan W inckel"2003' on post-AGB stars. 'Zuc kermanI 
on dusty disks and Waters & Waelkens 1998 on Herbig Ae/Be 
stars). 

However, the relatively large beam size of IRAS (45"x 15"at 
12/jm and larger at longer wavelengths) meant that regions such as 
the Galactic Plane could not be studied properly because of source 
confusion. In addition, due to its orbit, the IRAS satellite also did 
not observe 4% of the sky. These so-called "IRAS gaps" have never 
been surveyed in the mid-IR. 

The Mid-course Space Experiment (MSX) satellite carried 
out a mid-IR survey of the galactic plane and other areas of the 
sky missed by IR AS using the SPIRIT III instrument (full details 
IPrice et al.l200ll) . The sensitivity of SPIRIT III at 8^m is compa- 
rable to IRAS, but its beam size is approximately 35 times smaller 
(at the longest wavelength) and is therefore much less hampered 
by source confusion. It discovered 430 000 objects in the Galactic 
Plane (defined as | b |<6°), which is four times as many as IRAS 
detected in this region of the sky. Moreover, it surveyed the IRAS 
gaps for the first time in the mid-IR. 

The MSX survey thus provides an excellent opportunity to 
systematically search the Galactic Plane and the IRAS gaps for 
IR excess stars. The purpose of the present paper is to find opti- 
cally bright objects within the MSX Point Source Catalogues that 
were previously unknown as having IR excess. The resulting list 
can then be used for (optical) follow up studies. The methodol- 
ogy we adopt is similar to that of |Oudmaiier et al. 1 1992). They 
cross-correlated the optical SAO catalogue with the IRAS PSC, to 
identify objects with IR excess. The far-IR colours immediately re- 
vealed non-photospheric emission if the temperatures are less than 
about 200K while those objects with higher colour temperatures 
are identified by assessing the IR fluxes compared to that predicted 
from the photospheric optical emission. After a further selection 
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on spectral type (BAFG - typical for most evolved and young stars) 
this search revealed a master sample of 462 objects. As MSX did 
not observe beyond 20/im, it is not readily possible to find excess 
objects based on MSX data alone, and we will rely on computing 
the IR excesses using additional information to identify the stars of 
interest. Rather than use the SAO catalogue, we will use the much 
larger optical Tycho-2 catalogue ( H02 et al. 2000) and its accompa- 
nying spectral type catalogue li Wright et al.ll2003l) which contains 
more than 350 000 spectral types as a starting point. To further char- 
acterise the infrared emission of any objects found we will also use 
data from the 2MASS All Sky Catalogue of Point Sources (here- 
after 2MASS PSC) (Cutri et al. 2003). 

This paper is organised as follows. Section|2|describes the ini- 
tial input catalogues and the cross correlation of the Galactic Plane 
sample in detail. In Section|3|we discuss the properties of the re- 
sulting sample and develop a new method of identifying IR excess 
stars. We end with a discussion of the resulting final sample of ex- 
cess stars and outline the way it is presented. We also discuss in 
Appendix IXI the application of the excess identification procedure 
to the MSX IRAS gap survey and other regions surveyed by the 
MSX mission. 



2 INPUT DATA AND CROSS CORRELATION 

In order to define an IR excess for the sample stars, we use the MSX 
mid-IR catalogue, the optical Tycho-2 catalogue and the near-IR 
2MASS PSC. As will be shown later, once the "optical" stars have 
been identified, the use of their near-IR magnitudes has marked ad- 
vantages over using the optical data alone. Below we will describe 
the input catalogues and Galactic Plane cross-correlation in detail. 
The other MSX catalogue cross-correlations will be discussed in 
the appendix. 

2.1 Data Sources 

The MSX SPIRIT III instrument observed in six photometric bands 
between 4 and 21.34/im, and was most sensitive at 8.3/im where 
its sensitivity was comparable to the IRAS 12^m band. The MSX 
mission surveyed the galactic plane, IRAS gaps, Magellanic clouds 
as well as some other regions of high stellar density. A major data 
product of the MSX mission was the MSX PSC v2.3 (hereafter 
MSX PSC) I Egan et al. 2003) which has a positional accuracy, o — 
2", and a limiting flux at 8/im of 0.1 Jy. 

As we wish to identify a large sample of new IR excess stars, 
we require a sufficiently large and accurate optical star catalogue. 
For this reason we used the Tycho-2 star catalogue compiled from 
data collected by the Tycho star mapper on-board the Hipparcos 
astrometric satellite. The positional accuracy of Tycho-2 is excel- 
lent and more than required for our cross-correlation purposes. The 
completeness limit (99%) of the catalogue is given as 1 1 in the Vt 
band, and its limiting magnitude is about 14 in the Vt band and 15 
in the Bt band. Tycho-2 observed objects many times, resulting 
in a photometric precision of typically 0.05 mag. Obviously, this 
precision becomes worse at the fainter end, for magnitudes around 
10-11 the photometric error is 0.1 mag rising to 0.3 mag close to 
the detection limits. 

The Tycho-2 catalogue does not contain any stars brighter than 
Bt < 2.1 or Vt < 1-9 due to the nature of the data reduction 
technique. We used the brighter stars from the first supplement to 
the Tycho 2 stars catalogue. 

Near-IR photometry is taken from the 2MASS PSC which 



contains nearly half a billion sources. The survey is complete to 
magnitude K ~15. The data has relatively bright saturation limits 
(K ~3.5), though for A' < 8 the 2MASS catalogue lists photome- 
try estimated from radial profile fitting to the wings of the saturated 
sources. The photometric accuracy, especially at the fainter end, is 
better than that of Tycho-2 and for bright saturated sources com- 
parison with previous observations suggest the error is typically 
< 10%. 

Finally, to learn more abo ut the sample , spectral type infor- 
mation is very important. Wri ght et alJ j2003) compiled all known 
spectral classifications for Tycho-2 stars and their catalogue con- 
tains spectral types for 351 863 stars of which 61 472 are in the 
galactic plane (defined here as | b |<6°). As might be expected the 
objects with spectral types are somewhat brighter on average than 
the full sample. 

2.2 Positional Correlation 

We have positionally correlated the Tycho 2 star catalogue and the 
MSX v2.3 PSC. For a positional association between a Tycho 2 star 
and a MSX source, the separation of the two objects is required to 
be less than 6", corresponding to a typical 3a accuracy of the MSX 
position. In the case of two stars within 6" of an MSX source we 
have taken the closest. 

The deviations of the MSX positions compared to the Tycho-2 
positions in the right ascension and declination directions are cen- 
tred at zero with cr=2" indicating that, in a statistical sense, the 
Tycho-2 and MSX associations are real. The shape of the distribu- 
tion compared well with the similar in-scan and cross-scan separa- 
tion distributions shown by Egan et al. ( 2003|). To test the 6" cutoff 
the correlation was extended out to a radius of 15". We found that 
the majority (90%) of the associations had separations smaller than 
6" and the distribution showed that at larger radii the number of 
associations reached some random background level. We therefore 
concluded that 6" was the optimum value for the cutoff for posi- 
tional association. 

The MSX-Tycho 2 sources have also been further correlated 
with the 2MASS PSC, using a 6" search radius around the MSX 
position. In the case of multiple 2MASS sources within the 6" 
search radius, we have taken the closest match and flagged the 
MSX source. 

To ensure that all MSX-2MASS counterparts are also Tycho 
2-MSX-2MASS associations, we require the 2MASS source to be 
additionally within 0.50" of the Tycho 2 catalogue position. Failing 
to do the latter returns many sources that are either unrelated to the 
MSX or the Tycho-2 source and would severely contaminate the 
final sample. 

We made use of the Infrared Science Archive (IRSA)^ to per- 
form the 2MASS cross-correlations. 

2.3 Photometric Constraints 

At this stage we made certain requirements of the photometry pro- 
vided by the different catalogues for further inclusion in the sample 
so presence of the IR excess is credible/significant. For the Tycho 2 
star catalogue we required that a star be detected at both the Bt and 
Vt bands (corresponding to a magnitude less than 15 and 14 in Bt 
and Vt respectively according to the Tycho 2 explanatory supple- 
ment ,Hog_etaL„200^. For the MSX data we required a detection 

^ http://irsa.ipac.caltech.edu 
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12 783, are listed in the Tycho 2 spectral type catalogue. A smaller 
fraction of these 7443 have two dimensional spectral types. 



2.5 Comparison with IRAS 

To ascertain how many of our sources have previously been de- 
tected with IRAS and hence estimate how many objects in our fi- 
nal sample are new identifications, a positi onal correlation of the 
MSX-Tycho 2 sources with the IRAS PSC teei chman et alJl988^ 
has also been undertaken. We used the VizieR^ catalogue access 
tool to search the IRAS PSC around the MSX positions of our 
MSX-Tycho 2 associated stars. For the correlation we use a circular 
search radius of 45" typically corresponding to three times the ma- 
jor axis of the elliptically shaped IRAS io positional uncertainties 
(Beichman et al. 1988). 

Of the 35 044 MSX-Tycho 2 sources, 9473 were found to be 
within 45" of an IRAS source. We therefore conclude that the ma- 
jority of any IR excess stars in this paper are new identifications 
and previously unstudied. 



Figure 1. The magnitude distribution of the galactic plane sample in the op- 
tical and mid-IR. The optical By and bands are represented by the thick 
and faint solid lines respectively. The MSX 8/im distribution is indicated by 
the broken line. 



in at least one band with flux quality greater than 1 (signal-to-noise 
ratio greater than 5). The only constraint we made on the 2MASS 
photometry is that it must be better than an upper limit (i.e not hav- 
ing a U or X flag) at all bands. 



3 RESULTS FOR THE GALACTIC PLANE 

In the following we will first discuss the properties of the sample 
using colour-colour diagrams, and then continue with the identi- 
fication of excess stars. In order to do this, we will first derive a 
relationship between the near-IR colours and the MSX 8^m mag- 
nitude for normal stars. We will concentrate on the largest MSX 
catalogue, that of the Galactic Plane and briefly discuss the selec- 
tion from the other samples later. 



2.4 Galactic Plane 

The MSX PSC (ver 2.3) Galactic Plane | fe | ^ 6° survey (including 
the plane regions of the IRAS gaps) contains 43 1 7 1 1 sources. After 
cross correlation with the Tycho 2 optical star catalogue, we are 
left with 35 044 (« 8%) associations which meet our photometric 
constraints. The bulk (~ 75%) of these galactic plane sources were 
only detected at the most sensitive MSX band A (8.28^m). The 
further correlation of the MSX-Tycho 2 sources with the 2MASS 
PSC search found 33 495 (95%) counterparts. A small fraction of 
these (3%) had multiple 2MASS sources within the MSX search 
area, in these cases we cannot rule out the possibility that the MSX 
source is actually a nearby optically invisible 2MASS source rather 
than the Tycho 2 star. A clear example of this is the B type star 
HD 93942 (the object with K-[%] « 8 in Figure |4} which is in 
close proximity to a very red (optically invisible) carbon star. In the 
near-IR there are two 2MASS sources within 6" of the MSX source 
position. Although the MSX flux is most likely from the carbon 
star it is difficult to determine which of the 2MASS sources is the 
correct association. As we cannot objectively remove these sources 
we flag them. 

Using the magnitude distribution of the sample (see Figure^ 
we estimate the sample is complete to magnitude 8 in Vr, 1 1 in Bt 
and 5 in MSX Band A (or [8] hereafter^). This makes it the most 
complete IR sample of optical stars in the galactic plane to date. Of 
the 35 044 Tycho 2-MSX sources described above, about a third. 



The MSX fluxe s are converted to magnitudes using zeropoint fluxes from 
lEgan et alj l2003ft which is 58.49 Jy at 8^tm 



3.1 General Properties 

Figure |2| shows a (Bt — Vt,Vt — [8]) colour-colour diagram 
which is the MSX analogue to the well studied diag nostic IRAS 
{B-V,V-[12\) diagram (e.g. lWaters. Cote and Aumana.1987. her- 
after WCA). 

The left hand panel, containing all objects, shows a main band 
with increasing spread around it towards redder colours, the ob- 
jects with known spectral types are plotted in the right hand panel. 
Not surprisingly this subsample is on average brighter, and a large 
number of objects has dropped out, allowing us to recognise a well- 
defined band in the plot. The larger spread evident in the full sample 
is thus mostly due to larger photometric errors. 

This diagram is a good diagnostic to identify IR excess stars, 
as originally described by WCA. Normal stars follow a well- 
defined sequence, while objects with excess 8/im emission are 
readily identified by their deviating Vr — [8] colours. This is also 
observed here; the "main band" of stars follows a more or less well- 
defined relation and is accompanied by objects located above this 
relation. These are the IR excess stars, the majority of stars in the 
Tycho 2-MSX sample are not IR excess stars but normal stars. 

^ http://vizier.u-strasbg.fr 

On each diagram the magnitude and direction of a typical interstel- 
lar extinction vector are indicated calculated using the following assump- 
tions. For the mid-IR wavelengths (A > 5 ytm) we have used the MSX fil- 
ter averaged as t ronom ical silicate data lDraine & Leel 1 19841) as derived by 
iLumsden et'ai] |2002|) . For the near-IR wavelengths ( A < 5 ^im) we adopt 
an extinction law that varies as A"'^'^^. For the optical bands we have used 
the standard optical Ay = 3.1(-Eb-V') extinction relation. 
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Figure 2. Optical-mid IR colour diagram of the sources. Tlie horizontal axis denotes the Tycho Bt — Vt colours, while [8] represent the MSX 8/im magnitude 
(see text). The left hand panel shows all 35 030 MSX sources with a Tycho counterpart. The right hand panel shows the 12 778 objects with known spectral 
types. Note that this (brighter) sample shows a smaller spread around the main band. For comparison a reddening vector is shown as well. 
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Figure 3. Optical and Mid IR colour diagram showing the different regions occupied by different spectral type stars. The spectral type is shown in the lower 
right hand comer of the respective panel. The samples consist of 125 O-type stars, 1003 B stars. 718 A stars, 943 F stars, 2288 G stars, 6495 K stars and 1099 
M stars. Stars with General Catalogue of Galactic Carbon Stars identifications are also shown. The luminosity class of the stars is represented by the following 
symbols □ for I, II o III, IV and X for V, unclassified stars are indicated by a dot. The intrinsic colours of dwarf stars (see text for details) together with the 
extinction vector to indicate the degree of reddening of the sample are also shown in the upper left hand panel. 



3.2 Excess as function of spectral type 

To further investigate the properties of the sample, we plot the 
colour-colour diagrams for each spectral type in Figure |3| The 
sample is broken down into the spectral types 0-M while a set 



of Carbon st ars, taken from the Catalogue of Galactic Carbon 
Stars (CGCS lAlksnis et ai]|200l lh. is also shown. Where known, 
dwarfs (luminosity class V), giants (IV-III) and supergiants (I- 
II) are indicated by different plot symbols. As a guide for the 
intrinsic colours of main sequence stars in Figure |3| we show 
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the median Vr — [8] colours of BO to MO dwarf stars taken from 
ICohen. Hammerslev & Esm] i2000l) plotted against the intrinsic 
B—V colours of these stars given by Schmidt-Kalei ( 1982) which 
ar e corrected for Tycho photometry using the transformations given 
in lThe Hipparcos and Tvcho Catalogued 11997). 

The earlier type stars are dominated by main sequence stars, 
while giants constitute the majority for the later type objects as the 
dwarfs are too faint to be detected in both optical and IR catalogues. 
The transition between dwarfs and giants is visible as the relative 
paucity of sources at Bt — Vt ~ 0.75. 

A separate bifurcation is also visible at Bt~Vt ~ 1 in Fig|2| 
The upper sequence consists of reddened supergiants while the 
lower sequence consists of giant (and some dwarfs). This sepa- 
ration o f luminosity classes (beginning at early G type) was first 
noted bv lCohen et'ai]<1987h in the IRAS V- [12] diagram. A prop- 
erty apparent for all spectral types is that the non-excess stars span a 
much wider Bt ~ Vt range than the intrinsic colours for the respec- 
tive types signal. This is due to the presence of a large interstellar 
reddening at low galactic latitudes, and is illustrated by the red- 
dening vector indicated in the diagram. The main bands defined by 
non-excess stars per spectral type are parallel to the extinction, and 
correspond to values up to Av ~ 5. Often the reddest objects are 
supergiants, which is consistent with reddening as well; the objects 
with the largest extinction are presumably at the largest distances 
and therefore have to be intrinsically brighter. 

For the B type stars, and also, but less obvious, for the O stars 
due to the lower number of stars involved, the colour diagram is 
dominated by two bands which are both populated by Main Se- 
quence stars. The lower band contains normal, non-excess objects, 
and the upper band is populated by a large number of excess stars, 
many of which are Be stars. The excess emission is explained by 
free-free emission from the io nised gas in the circumstellar disk 
iWaters. Cote & LamerJl987l) . The relative number of A type ex- 
cess stars is much smaller, and seems to consist of dwarfs and su- 
pergiants in roughly equal numbers. There would seem to be a spe- 
cial class of A type stars around Bt~Vt = 1.5, where these seem to 
branch off at a different slope than the other A type stars. It appears 
however that these objects have colours similar to K type stars and 
the most likely explanation for this is that their spectral types are 
misclassified or that they are binaries with a K type counterpart. 

The F, G and K type stars have fewer excess objects. The M 
star sample shows a distinct upturn. As already noted by WCA, for 
these cool objects, the optical colours probe the Wien part of the 
spectral energy distribution and show hardly any dependency on 
temperature anymore. The B/im band is sensitive to spectral fea- 
tures, such as molecular bands and therefore shows a strong gradi- 
ent. This particular colour diagram is thus not a useful diagnostic 
to identify M stars with IR excess. 

Perhaps unsurprisingly, the carbon stars follow a similar band 
as the M stars. However the relatively large photometric uncertain- 
ties of these faint stars means they sometimes appear at much bluer 
colours than would normally be expected. 

It is important to highlight a significant dif ference with pre- 
vious studies here. WCA in-particular, but also lOudmaiier et alj 
il992h used optical catalogues w ith brighter cut-offs (the Bright 
Star Catalogue Hof fleit et ail99l | and the Smiths onian Astrophys- 
ical Observatory star catalogue's AO Staff'1966' with limits V = 7 
and 9-10 respectively) than used here. The deeper Tycho-2 sam- 
ple inevitably contains a large fraction of (heavily) reddened stars, 
which is aggravated by the fact that we observe in the direction of 
the Galactic Plane. It can be seen in the upper left panel of Fig- 
ure|3| that the extinction vector is non-parallel to the "normal star" 



relation, this gives heavily reddened normal stars Vt — [8] colours 
similar to excess stars. Moreover, the fainter, numerous, Tycho-2 
stars have comparatively large photometric errors (as indicated by 
the large spread in Figure|2}, making it harder to recognise excess 
emission from the optical-IR colour-colour diagram. It is therefore 
not trivial to properly identify stars with excess from the Tycho 2- 
MSX data alone. Therefore, to make headway, we need to go to 
wavelengths where the extinction will be less, the near-IR. 

3.3 Near-infrared 

To reduce the impact of extinction on the IR excess star selection, 
we show the near- and mid-IR colour equivalent to Fig|2|in Fig.|4] 
The K— [8] range spanned by the normal stars is much smaller than 
the Vt — [8] range because we now probe the Rayleigh- Jeans tail 
of the spectral energy distribution for all stars. The H—K colour 
range is small for the same reason, while the reddening is less se- 
vere (typically E{H-K) = 0.22E(B~V)). Most of the objects 
beyond H —K>Q.5 suffer from reddening, also objects with excess 
emission at K give rise to red H—K colours. 

To understand and describe the features of the colour di- 
agram, we have also plotted the colours of known objects 
(FiglU which have been observed by MSX (but are not nec- 
essarily part of the Tycho 2 sample). The objects plotted are 
from the following c atalogues: Optical selected carbon stars from 
lAlksnis et 51 <200lh , OH/IR stars (AGB stars with large mass 
loss rates stars) from [ Chen g^ur et a i iT995) . Herbig Ae/Be stars 
from lThe. de Winter & Perea <I994 and the Be type stars extracted 
from our own Tycho 2 sample (which were found to agree with 
the mid- near- IR colours of all known Be stars as compiled by 
IZhang. Chen'~& Y ang 2005). 

The distribution of B stars with excess in Fig.|5]shows a clearly 
defined band of excess stars with if— [8] < 2 and a more diffuse 
distribution of stars with greater excess emission, which is further 
complemented by a large number of excess stars without spectral 
classification (see Fig.|4j. 

It is immediately apparent from Fig|4|that only the B stars 
with -R"— [8] < 2 can be explained by the properties of the known 
Be stars. The remaining objects with greater excess emission are 
therefore very interesting; we discuss in section|5|the possible na- 
ture of these objects. 

The two branches beginning H—K = 0.50, with gradients 
roughly equal to the reddening vector , are clearly associated with 
the carbon and OH/IR stars. The OH/IR stars also have a slightly 
greater (1 mag) -ft"— [8] colour, presumably due to their higher mass 
loss rates over these optically selected carbon stars. The Herbig 
Ae/Be stars all show excess emission and are heavily reddened, 
due to their surrounding circumstellar material. 

Spot checks on stars with H—K < revealed the majority 
have saturated K band magnitudes, and therefore their colours are 
just affected by the resulting errorbar. The objects with -?f— [8] < 
—0.50 are in general variable stars and binaries. 

Although we plot the same number of objects in Figs|2|and|4] 
the excess stars are more easily distinguishable than in the optical 
diagram. In addition to the reduced extinction, this is also the case 
because of the smaller photometric errorbars involved, validating 
the use of near-IR colours. 

The distribution over spectral type is shown in Figure |5] All 
features described for the optical (Fig|3} are present, but much more 
prominent. It is clear that the near-mid IR two colour diagram is a 
much better tool for the selection of IR excess stars, due to the 
removal of the strong extinction effect and the better photometric 
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Figure 4. Near-Mid IR (left) colour diagram of MSX-Tycho2-2MASS sources (33 485), with spectral type (centre) (12 452) and (right) the colours of known 
sources (not necessarily in our sample). The excess cutoff line is shown in the left panel (see text for details). The symbols in the right hand panel indicate the 
type of objects (* Be type star from our sample, x a OH/IR star.n a Herbig Ae/Be star and a dot indicates an optical Carbon Star). 




Figure 5. Near-Mid IR colour diagram showing the different regions occupied by different spectral type stars. The spectral type is shown in the upper left 
hand comer. The samples consist of 120 0-type stars, 937 B stars, 672 A stars, 885 F stars, 2219 G stars, 6406 K stars and 1086 M stars. Stars with CGCS 
identifications are also shown.The luminosity classification of the stars is represented by the following symbols □ for o III,IV and x for V, unclassified 
stars are indicated by a dot. The intrinsic colours of dwarf main sequence stars (see text for details) together with the extinction vector to indicate degree of 
reddening of the sample are are also shown in the lower right panel. 
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accuracy of 2MASS for the optically fainter objects. We therefore 
proceed with this colour diagram to identify IR excess stars. 

3.4 The H-'K, K—[%\ relation for normal stars 

To identify stars with excess 8/im emission, we first need to know 
the photospheric contribution for all stars at this wavelength. To 
derive this we adopt the iterative procedure used by WCA for IRAS 
sources. This involves dividing the sample into H—K colour bins 
(0.01 mag width in our case) and calculating the mean {K~\S\)av 
and standard deviation a of the bin. By eliminating the outliers 
(defined as deviating by more than 2a from the mean), this process 
was repeated until the iK~[8\)av no longer changed significantly 
with further iterations. 

The procedure was only ru n over the intrins ic main se- 
quence near-IR colours given by iKoornneel il983 ') (H—K — 
—0.05.. 0.45). This step essentially removes the heavily extincted 
sources, whose infrared colours are heavily distorted by reddening. 
The final {K— [&\)av values could then be well fitted by the follow- 
ing straight line. 

{K ~ m^hoto = QAl{H - K), + 0.05 (1) 

This photospheric relation was extrapolated to cover the en- 
tire range of {H~K) colours spanned by the sample stars. Using 
thi s relation an d the intrinsic main sequence near IR colours quoted 
bv lKoomneed Jl983l) we derive a relationship between the spectral 
type of a star and its K— [8] colour (as shown in the lower right 
panel of Figure |5j. 



4 TYCHO-2 STARS WITH INFRARED EXCESS 

In this section we select the excess stars. We first discuss the more 
numerous 8/im excess stars, and then discuss excess at the other 
MSX wavelengths. 

4.1 Stars with warm dust : Excess 8 ^m emission objects 

To calculate the IR excess emission of the stars in our sample we 
follow O udmai jer et al. ( 1992) who used the relationship between 
the photospheric colours in the optical and IRAS 12 fj,m for normal 
stars. With the expression of the intrinsic ii"— [8] colour for nor- 
mal stars already in hand, it is trivial to compute the excess 8^m 
emission. Expressed in magnitudes this is done by E{K—[8]) — 
(A'-[8])-(if-[8])photo. 

The main issue is to decide on where to choose the cutoff for 
inclusion into the final sample of infrared excess stars. In Fig. |6| 
we show a histogram of the excesses in magnitudes. The distribu- 
tion of E{K~[8]) is asymmetric about zero. The negative excesses 
show a Gaussian shape. The positive side also displays an underly- 
ing Gaussian distribution, with a large non-gaussian tail that con- 
tains the excess stars. The gaussian shape on the negative side is 
explained simply by the width of the observed sample in the H—K 
A'- [8] colour diagram. The width (FWHM = 0.20) in this case 
is dominated by the photometric error in the MSX 8/im flux (5- 
20%) rather than the intrinsic scatter of normal non-excess stars. 

We can now derive the fraction of objects with excess 8^m 
emission. By using the negative excess distribution (assumed to be 
symmetric around zero) as an estimate of the non-excess distribu- 
tion, we can calculate the percentage of excess stars as a function 
of E{K—\8\). This is shown in Figure|6| the percentage of excess 
stars is zero for the negative side (by definition) and rises relatively 



quickly to 100% at E{K-[8\) = 0.80. We choose to have at least 
an 80% probability of IR excess in our sample. This corresponds to 
a E{K-[&) cutoff of 0.40 magnitudes. 

Using this criteria we identify 1830 8/im excess stars in our 
sample, around half of which (965) are not within 45"of an IRAS 
Point Source and are therefore new identifications. We also identify 
one excess star which was too bright to be included in the Tycho 
2 catalogue (and hence was picked up in the Tycho 2 supplement), 
the well known Be star 7 Gas. 

4.2 Stars with hot dust : K band excess emission objects 

It should be noted that the above procedure selects objects that have 
excess 8/im emission relative to the K band. An object surrounded 
by only hot dust will not be selected as its 8/im emission has a sim- 
ilar A— [8] colour as a cool C or M type star. Indeed, the intrinsic 
A- [8] colour for a 1000-2000 K Black Body is about ~ 0.5, and 
it may well be that such objects will not be selected as their ob- 
served Jf— [8] colour seems to all intents and purposes "normal". 
We therefore need to identify K band excess stars separately. 

We investigated several ways to identify such objects and 
found the best method is to use Vt — J, J— [8] colours as shown 
in Figure0 This avoids the problem of having the large errorbars 
we encountered in the optical Bt — Vr colours and includes objects 
with excess even at the H band. It is extremely rare for hot dust to be 
responsible for J band excess. Indeed, checks on a random sample 
of / excess emission objects revealed them to be uncatalogued bi- 
nary objects, where the 7 emission is due to a bright cool secondary. 
We identify the K band excess stars using a similar iterative proce- 
dure as for the A— [8] colours and we derive the following relation 
for the photospheric contribution. 

(J - mphoto = -0.025 + QA2Q{Vt - J) - 0.021(Vt - J)\2) 

The excess J— [8] colour distribution showed that the proba- 
bility of excess reached 80% at a E{J—\8]) = 0.60 magnitudes 
and this is where we placed the excess cutoff. This method of ex- 
cess identification returned only 68% of the A— [8] excess sources 
and a further 95 stars which do not have excess compared to K at 
8/im or 21/im. Interestingly only 9 were not observed by IRAS, 
possibly because the hot excess stars are relatively infrared bright. 

4.3 Stars with cool dust : 21 ^m excess emission objects 

Stars such as post-AGB stars surrounded by a cool detached dust 
shell may not show any distinguishing excess emission at 8/jm. We 
therefore wish to use the longer wavelength MSX bands to identify 
such source s in our sample. The pr evious searches for cold excess 
objects (e.e lOudmaiier et al]|l992h made use of the IRAS far-IR 
colour-colour diagram to identify excess if temperatures are less 
than 200K (indicating that infrared emission was not photospheric). 
A similar approach with MSX is hampered by the reduced sensitiv- 
ity of MSX at longer wavelengths. Indeed, very few sources (2%) 
in the sample are detected (i.e flux quality greater than 1) in every 
band. As a compromise we will use the K— [21] colour (see Figure 
to identify cool excess objects. This has the further advantage 
that stars detected exclusively at 21^m are not removed from the 
sample as they would be if we relied on a multi MSX colour dia- 
gram. 

FigureQshows the A— [21] colours, of all sources with a flux 
quality greater than 1 at 21/im, plotted against H~K. Objects that 
are found to display 8^m excess emission are indicated by a dot. 
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Figure 6. Excess K — [8] {E{K — [8]) distribution of sample, tlie negative excess distribution is mirrored tlirougli E{K — [8]) = to illustrate the asymmetry 
of the distribution. The right hand panel shows the percentage of excess stars as a function of E{K — [8]) emission, the excess definition threshold is also 
shown at 0.40 magnitudes. 







^ 1 1 1 1 1 1 I 1 1 1 


1 1 1 1 1 1 1 1 1 1 4 




B^llPinO- - 


















1- 1 1 1 1 1 1 1 1 1 





-2 2 4 6 
V-J 



10 



10 



5 - 



- 




5A„ 



0.5 1 

H-K 



1.5 



Figure 7. Optical Mid infrared and Near-Mid infrared colour diagrams showing the J- [8] and K-[21] colours for hot and cold excess emission respectively. In 
the left panel, a □ indicates a star with excess J-[8] emission but normal K-[8] colours, the rest of the sample is indicated by a dot. The right panel shows the 
K-[21] colours for cool excess identitication. Stars detected at 21^m but not at 8/xm are indicated by a *, stars detected at both 8fim and 21/im are indicated 
by a dot if they show excess at 8/xm or a X if they do not 



while those not detected at 8^m or not found to have excess are 
indicated by larger plot symbols. It is immediately apparent that the 
majority of the cool excess objects also have excess %[im emission. 

Of the stars without 8/im excess only a few have if— [21] 
colours significantly different from the rest of the sample (the ma- 
jority of which are M giants). All sources that are detected at 
21/im but are not detected in the more sensitive 8pm band show 
cool excess emission. To determine the ex cess cutoff we have used 
the median V— [25] values published by ICohen et alj il987h and 
converted t o A"— [25] using the appropriate V—K colours from 
iKoomneed h983h . As the [21] and [25] bands both probe the 
molecular absorption bands of the cool M type giant stars, we may 
assume that the if— [25] and if— [21] colours are roughly com- 
parable. We therefore define the excess cutoff to be the 3 sigma 



deviation from the intrinsic if— [25] colours of giants. Hence we 
select stars with cool infrared excess if they have H—K < 0.1 
and K-[21] > 1 or if- [21] > 3 at H-K > 0.1, as shown in 
Fig0 Using this criteria we identify 13 cool IR excess stars which 
do not have excess 8^m emission, 5 of which do not have IRAS 
counterparts. 



5 DISCUSSION 

In this paper we have searched for stars with infrared excess in the 
combined Tycho-2, 2MASS and MSX catalogues. This resulted in 
a grand total of 1938 excess stars, of which 979 were not detected 
by IRAS and are thus newly discovered objects. 
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5.1 Distribution over Spectral Type 

In Figure|8|we show the number of Tycho 2-MSX stars as a func- 
tion of temperature class. Alongside we also show the percentage 
of these stars which were found to have excess emission. 

We find that the percentage of excess appears to decrease with 
effective temperature and that the hotter stars (0,B and A) and very 
cool (M) have a much higher percentage of excess stars than the F,G 
and K classes. However it is interesting to note that the numbers of 
Tycho 2-MSX stars detected shows the exact opposite behaviour, 
indicating that a strong selection effect may be present. 

A comparison of the Tycho 2-MSX sample with the entire 
galactic plane Tycho 2 spectral type catalogue shows that we de- 
tect a higher fraction of G and K type compared to other spectral 
types. The most obvious source of this effect is the sensitivity of 
the MSX satellite. For example, the photospheric emission at 8/im 
for all but the brightest 0,B and A type stars may be below the 
detection threshold of MSX; we only detect the excess and bright 
non-excess stars (which are relatively few) thus leading to a higher 
fraction of excess stars for these spectral types. 



5.2 B type stars with infrared excess 

The high fraction (30%) of excess in B type stars is surprising, as 
the fraction of Be type stars in the B right Star Catalogue is well 
known to be 15% ICote and van Kerkwiik 1993) . We also noted 
earlier that there exists a large group of B type (with and without 
spectral classification) excess stars with colours that were signifi- 
cantly different to the known Be stars. The nature of these objects 
is not altogether obvious. 

Their H—K colours are quite blue indicating that these are in- 
deed B type stars. The normality of the H—K colour indicates that 
they do not have a great deal of K band excess. The absence of ex- 
cess in the K band would seem to lead us away from a free-free 
emission explanation for these excesses and indicates that these 
stars are unrelated to the Be stars immediately below them in Figure 
|4| The objects are also not particularly reddened, ruling out heavily 
dust embedded objects and implying the surrounding dust must be 
relatively optically thin. 

We are therefore looking for a optically thin thermally ra- 
diating dust mechanism for these stars. Possibilities include: Hot 
Post-AGB stars with a optically thin cool detached dust shell, re- 
flection nebulosities, giant stars heating the surrounding interstellar 
medium dust, weak HII regions and Vega type stars surrounded by 
a close warm disk of dust. 

To resolve the situation we looked at the spectral energy dis- 
tributions of the subsample with IRAS detections. We showed that 
the objects with intermediate -ft^— [8] colours (2 — 4) to be Herbig 
stars. Be stars, reflection nebulae and post-AGB stars. The objects 
with larger K— [8] colours were found to be predominately objects 
such as post-AGB stars, HII regions and a number of galactic cirrus 
contaminated fields (some possibly also heated by giant stars). We 
did not find much evidence for Vega type disk systems. To com- 
plement this we also inspected Galactic Legacy Infrared Mid-P lane 
Survey Extraordinaire (GLIMPSE; see Benja min et im- 
ages for the small number (~ 10) of these stars which were within 
the survey region. The GLIMPSE images mostly showed these ob- 
jects to be stars within diffuse emission regions and/or associated 
with star forming regions. These objects would benefit from ground 
based optical follow up observations to better confirm their nature 
and to further understand this region of the colour diagram. 



5.3 Presentation of the Data: Data Tables 

We publish in hardcopy (see Table0 the stars which were not de- 
tected by IRAS (and therefore should be new identifications), that 
have full spectral type information and excess — [8] > 0.75 mag 
and/or excess emission at other bands. The entire sample (includ- 
ing the other regions discussed in the Appendix) is published in the 
online version of this paper (see Tables 2,3,4 and 5). A machine 
readable format will also be made available at CDS Strasbourg . 

5.4 Final Remarks 

To identify stars with IR excess emission we have positionally cor- 
related the MSX PSC with the Tycho 2 star catalogue. We found 
that a near-mid IR colour diagram had marked advantages over a 
optical-mid IR colour diagram in selecting these objects, specif- 
ically in reducing the strong line of sight extinction experienced 
in the direction of the galactic plane. Using the derived colours of 
normal stars a selection criteria was developed to identify excess 
emission sources from the colour diagram. The criteria produced a 
sample of 1938 stars in the galactic plane, just over 50% of which 
were determined to be new identifications of infrared excess (i.e 
not previously detected by IRAS). 

The majority of the excess stars were found to be hot stars, a 
high (30%) fraction of B type stars with excess was discovered. The 
infrared excesses for a number of these B type stars were found to 
be much greater than those of the known Be sample. The known ob- 
jects in this group were found to be a mixture of Herbig stars, post- 
AGB stars, reflection nebulae and stars contaminated by galactic 
cirrus. 

The other regions surveyed by MSX were also searched for 
IR excess stars and the same selection criteria were applied (see 
Sec|^. We publish the entire excess sample online (see Table 2) 
and the brightest new identifications of excess stars in the galactic 
plane, in Table Q These lists should provide a good starting point 
for future ground based follow up observations. 
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Figure 8. Temperature class distribution of entire MSX-Tycho 2 sample (left) and percentage of stars with excess emission as a function of spectral type 
(right). 



Notes to Table 1 

Table 1 contains: the Tycho identifier, Tycho (J2000) position, 
HP name f rom the HD catalogue identifications for Tycho 2 stars 
iFabriciu~et aL.2002'). sp ectral type from the Tycho 2 spectral type 
catalogue <Wrightetal1 l2003) , Tycho Bt and Vt band magni- 
tudes, 2MASS J,H,K near-IR magnitudes, a flag indicating the 
number of extra 2MASS sources within 6" of the MSX position 
(# 2MASS), the MSX 8,12,14 and 21 band fluxes (in Jansky), the 
photometric quality of each of the MSX bands (upper limit=l to ex- 
cellent=4) and finally the excess colour in magnitudes. The excess 
colour listed is -R"— [8] unless the entry has superscript J or 21, in 
which case it is excess J— [8] or if— [21] respectively. The typical 
errors in the Tycho 2 photometry are 0.05 mag rising to 0.3 mag 
close to the detection limits. The 2MASS photometry is of higher 
accuracy than Tycho 2 especially at the fainter end. 



Notes to the Online only Tables 2,3,4,5 and 6 

The tables 2,3,4,5 and 6 are only available in the online version of 
this paper. They have the same format as Table 1, aside from an 
extra column indicating the presence of an IRAS counterpart. The 
entry contains an T if an IRAS source was found within 45"and is 
empty if none was found. Additionally Table 6 also includes a col- 
umn indicating which star forming region the source is associated 
with. 



Table 1: Tycho 2 stars with infrared excess 



TYCl 


TYC2 


TYC3 


Ra 


Dec 


HD 


SpType 


Bt 


Vt 


J 


H 


K 


# 


F8 


F12 


F14 


F21 


Flux 


Excess 








(J2000) 


(J2000) 






mag. 


mag. 


mag. 


mag. 


mag. 


2MASS 


Jy 


Jy 


Jy 


Jy 


Quality 


mag. 


3664 


192 


1 


00 11 37.1 


+58 12 42 
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B5I1: 


7.26 


7.11 


6.49 


6.43 


6.33 





0.39 








4000 


0.79 


4034 


1319 


1 


01 24 19.4 


+62 49 36 




B1:PE(V) 


11.51 
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8.13 
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1 


0.15 








3 000 


1.17 


4032 


17 


1 


01 58 39.1 


+60 37 43 




B 1 VE 


10.81 


10.23 


8.32 
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7.74 





0.16 








3 


1.15 


3697 


1551 


1 


02 02 36.4 


+59 41 17 


12302 
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8.37 


8.12 


7.22 


7.09 


6.90 





0.33 








4 


1.16 


4041 


1668 


1 
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+64 01 28 


13590 
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8.33 


8.00 


7.14 


6.99 


6.83 





0.29 








3 


0.93 


3694 


1707 


1 


02 22 06.4 


+57 05 25 




B 2 ni-IVE 


10.05 


9.83 


8.48 


8.10 


7.55 





0.29 








4000 


1.53 


4049 


1502 


2 


03 16 17.4 


+60 02 07 


20053 


B 3 V 


8.43 


8.31 


7.98 


8.02 


8.00 


1 


0.17 


2.96 






2 100 


1.61 


3718 


756 


1 


04 04 21.6 


+53 19 44 


25348 


BIVPNNE 


8.44 


8.20 


7.34 


7.16 


6.94 





0.23 








4 


0.77 


2401 


381 


1 


05 09 56.4 


+37 00 16 


33152 


B 1 VE 


8.37 


8.15 


7.69 


7.66 


7.53 





0.25 








3 


1.51 


2900 


503 


1 


05 10 48.2 


+41 00 10 


33232 


B 2 VNE 


8.29 


8.25 


7.91 


7.81 


7.64 


1 


0.14 








2 


0.97 


2900 
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1 


05 13 13.3 


+40 1 1 36 


33604 


B2V:PE 


7.37 


7.39 


7.44 


7.46 


7.49 





0.18 








3 000 
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16 


1 
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1 
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7.95 


1 
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4.57 
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19.67 


4444 


3.89 
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1417 


1 
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+26 26 43 


39340 
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8.18 
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7.26 





0.23 
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7.62 


7.70 
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7.71 





0.11 
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0.81 


5381 


2964 


1 
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7.96 


7.76 


7.13 


7.10 


6.92 





0.23 
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0.79 
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3009 


1 


07 11 20.8 
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7.21 


7.31 
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1 
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1 
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u 


U.ZZ 








Q n n n 
U U U 


U. /D 


8645 


2059 


1 


12 14 01.8 


-59 23 48 


106309 


B2 111/VNE 


7.78 


7.86 


7.82 


7.83 


7.62 





0.18 








2000 


1.17 


8994 


1920 


1 


13 10 20.5 


-62 41 18 


114214 


FO V 


10.15 


9.39 


7.73 


7.50 


7.35 





0.19 








3 000 


1.00 


8994 


2766 


1 


13 14 11.5 


-63 22 25 


114800 


B2 ni/VNE 


8.05 


7.98 


7.51 


7.36 


7.10 





0.25 








4000 


1.01 


8994 


3012 


1 


13 20 35.5 


-63 24 43 


115746 


B2/4ni/V 


9.62 


9.38 


8.08 


7.92 


7.67 





0.14 








2000 


0.95 


8990 


3680 


1 


13 22 37.9 


-60 59 18 


116087 


B3 V 


4.34 


4.49 


5.24 


4.90 


4.89 





0.61 








4000 


0.65-' 



Table 1: (continued) 



TYCl 


TYC2 TYC3 


Ra 


Dec 


HD 


SpType 


Bt 


Vt 


J 


H 


K 


# 


F8 


F12 


F14 


F21 


Flux 


Excess 






/ r\f\r\\ 

(J2000) 


(J2000) 






mag. 


mag. 


mag. 


mag. 


mag. 


2MASS 


Jy 


Jy 


Jy 


Jy 


Quality 


mag. 


8995 


3220 1 


13 27 25.1 


-62 38 56 


116781 


09/Bl (I)E 


7.74 


7.67 


7.22 


6.99 


6.72 





0.44 








4000 


1.23 


9003 


1626 1 


13 28 07.3 


-66 16 47 


116849 


Bl VPE 


9.33 


9.17 


8.72 


8.52 


8.29 





0.13 








2000 


1.49 


8999 


53 1 


13 29 54.6 


-65 30 07 


117111 


B2 (V)NE 


7.78 


7.72 


7.28 




7.04 





0.24 




0.51 




4 10 


0.98 


9016 


519 1 


13 45 18.4 


-66 45 16 


1 19423 


B3/5 VNE 


7.54 


7.59 


7.59 


7.62 


7.59 





0.12 








2000 


0.81 


8676 


1771 1 


13 50 26.1 


-59 44 52 


120330 


B2/3 VNNE 


7.93 


7.91 


7.51 


7.41 


7.22 





0.25 








4 


1.16 


9005 


873 1 


14 15 09.9 


-61 06 42 


124340 


AO V 


8.42 


8.26 


7.77 


7.75 


7.68 





0.14 








2000 


1.02 


9007 


747 1 


14 35 33.3 


-61 00 29 


127756 


Bl/2 VNE 


7.69 


7.60 


7.06 


6.93 


6.66 





0.34 








4000 


0.92 


8692 


1383 1 


14 42 05.6 


-58 21 34 


128937 


B3/5 (III) 


7.88 


7.72 


7.34 


7.34 


7.24 





0.35 








4000 


1.58 


9024 


980 1 


14 59 21.9 


-62 30 32 


131925 


A3 IV 


8.51 


8.20 


7.45 


7.34 


7.23 





0.19 








3 000 


0.89 


9025 


1827 1 


15 09 03.3 


-61 53 15 


133738 


Bl/2 III/ VNE 


6.93 


6.97 


6.66 


6.62 


6.41 





0.47 








4 


1.04 


8702 


636 1 


15 10 59.6 


-57 42 43 




9.5IA 


11.51 


10.54 


7.65 


7.35 


7.12 





0.22 








4000 


0.90 


8706 


14 1 


15 15 16.2 


-58 10 22 


134958 


B2 (II)NE 


8.48 


8.20 


6.76 


6.57 


6.34 





0.44 








4000 


0.88 


8695 


2105 1 


15 25 55.3 


-53 46 15 


136968 


B5 VNE 


8.26 


8.25 


8.34 


8.28 


8.14 


1 


0.12 


1.85 






2 10 


1.33 


8304 


641 1 


15 42 19.8 


-49 59 52 


139790 


B2 (III)N 


8.62 


8.57 


7.92 


1.16 


7.55 





0.16 








2 


1.02 


8321 


711 1 


15 46 47.4 


-52 08 46 


140605 


B5/7 IV 


7.01 


7.07 


6.96 


6.95 


6.86 





0.34 








4 


1.19 


8697 


718 1 


15 55 29.6 


-52 48 15 


142170 


B8/9 (IV) 


9.95 


9.92 


9.57 


9.05 


8.54 


3 


0.13 








2000 


1.62 


8322 


952 1 


16 04 36.2 


-51 38 36 


143751 


B3 n/in 


10.23 


9.72 


8.56 


8.45 


8.36 





0.28 








4000 


2.48 


8323 


688 1 


16 12 43.0 


-51 35 54 


145300 


FO III 


10.91 


10.44 


9.46 


9.29 


9.21 





0.21 








4 


3.00 


8319 


698 1 


16 18 39.4 


-49 24 49 


146444 


B2 VNE 


7.67 


7.60 


7.15 


7.09 


6.93 





0.29 








4000 


1.04 


8715 


1647 1 


16 19 14.2 


-54 57 42 


146463 


B3 VNNE 


8.05 


8.08 


7.92 


7.90 


7.71 





0.15 








3 000 


1.11 


8324 


1544 1 


16 19 42.8 


-51 02 04 


146630 


K2 III 


8.57 


6.76 


3.96 


3.21 


3.06 





3.53 


2.11 


3.76 


12.00 


4444 


0.39^^ 


8333 


1004 1 


16 34 43.5 


-49 33 09 


330950 


BIVNE 


10.05 


9.78 


8.66 


8.56 


8.50 





0.14 




1.01 




20 10 


1.85 


7862 


115 1 


16 35 48.7 


-42 07 22 


149313 


B2 lA/BE 


9.95 


9.48 


7.95 


7.72 


7.48 





0.17 








3 


0.97 


8333 


1429 1 


16 36 11.8 


-49 15 47 


149298 


B2 (II) 


10.25 


9.77 


8.12 


7.91 


7.65 





0.16 








2000 


1.06 


8329 


1025 1 


16 38 57.7 


-47 24 02 


149742 


K2ni 


10.23 


8.72 


6.22 


5.60 


5.43 





0.99 


0.94 


0.48 




4110 


0.87 


8338 


2080 1 


16 47 51.4 


-51 46 04 


151083 


B2 VN 


9.25 


9.08 


8.21 


8.04 


7.76 





0.15 








2000 


1.15 


7868 


448 1 


16 51 00.0 


-37 30 52 


151771 


B8 n/in 


6.33 


6.23 


7.79 


7.72 


7.66 


1 


0.32 








4000 


1.92 


8335 


1898 1 


17 00 28.7 


-49 15 14 


153222 


Bl IB/IIIE 


9.17 


8.90 


8.38 


8.30 


8.18 





0.15 








2 


1.61 


7372 


375 1 


17 05 52.8 


-36 35 17 


154243 


B3 VNNE 


8.28 


8.08 


7.29 


7.14 


6.92 





0.29 






2.42 


400 1 


1.02 


7878 


1149 1 


17 16 17.6 


-42 20 20 


155896 


B2/3 (V)NNE 


7.03 


6.96 


6.20 


6.12 


5.99 





0.63 








4000 


0.97 


7374 


641 1 


17 18 38.3 


-36 05 13 


156369 


A2 V 


8.94 


8.32 


6.76 


6.61 


6.44 





0.37 








4000 


0.83 


7874 


809 1 


17 19 16.8 


-39 48 25 


156409 


B2 (II)NE 


9.18 


8.78 


7.94 


7.80 


7.59 





0.18 








3 


1.20 


7384 


16 1 


17 37 30.9 


-35 19 59 


159684 


B2 VNE 


8.63 


8.24 


6.80 


6.60 


6.34 





0.53 








4 


1.07 


6840 


1414 1 


17 48 35.5 


-29 57 28 


316341 


B0,5V(PE)? 


9.74 


9.15 


7.60 


7.34 


6.96 


1 


0.26 








4000 


0.85 


7377 


827 1 


17 49 08.6 


-31 17 16 


161839 


B5/7 Win 


9.77 


9.64 


7.63 


6.58 


6.35 


1 


1.33 


0.89 


1.04 




4 12 


2.09 


6841 


1403 1 


17 59 47.6 


-23 48 58 


163955 


B9 V 


4.70 


4.72 


4.84 


4.61 


4.47 





1.29 


0.72 


0.68 




4 120 


0.78'^ 


1 JOZ 


1^1/1 1 

1 J 14 1 


1 1 '^Q A 
L / jy JO. 4 




10.:>(5Dc5 


DZ/j { V )I\tl 


/ .jV 


/.JJ 


1 . 14 


/.UD 


D.oo 


n 

u 


n ^ 1 
\}.j 1 








^ A A n 
J U U U 


1 .u / 


6846 


437 1 


18 04 39.6 


-26 01 07 


164950 


B3/5 II 


9.43 


9.13 


8.07 


8.00 


7.77 





0.13 








2000 


1.01 


6263 


2940 1 


18 04 43.2 


-20 56 44 


165014 


F2 V 


10.11 


9.33 


7.42 


7.04 


6.77 





1.08 


0.78 


0.76 




4 110 


2.27 


6259 


2481 1 


18 05 58.8 


-19 5713 


165285 


Bl/2 (I)NN(E) 


8.89 


8.45 


7.02 


6.78 


6.48 





0.44 








4000 


1.00 


6263 


2328 1 


18 07 11.4 


-21 26 38 


165516 


Bl/2 IB 


6.32 


6.26 


6.00 


5.99 


5.90 





0.60 


0.88 






4 100 


0.83 


6272 


2074 1 


18 08 27.1 


-19 52 07 


165783 


B3/5 IB 


8.69 


8.37 


7.45 


7.28 


7.09 





0.22 








4000 


0.92 



Table 1: (continued) 



TYCl 


TYC2 TYC3 


Ra 


Dec 


HD 


SpType 


Bt 


Vt 


J 


H 


K 


# 


F8 


F12 


F14 


F21 


Flux 


Excess 






(J2000) 


(J2000) 






mag. 


mag. 


mag. 


mag. 


mag. 


2MASS 


Jy 


Jy 


Jy 


Jy 


Quality 


mag. 


6268 


2490 1 


18 10 18.3 


-18 11 41 


166188 


B3 (II)E 


9.31 


8.98 


7.68 


7.52 


7.28 





0.23 








4000 


1.09 


6268 


212 1 


18 10 38.8 


-17 44 02 


166288 


B5 IB 


10.38 


10.06 


9.06 


8.99 


8.92 





0.28 








4000 


3.03 


6268 


620 1 


18 10 42.1 


-17 55 05 


166289 


Ki/2 (in) 


11.17 


9.31 


6.35 


5.64 


5.45 





1.45 


1.11 


1.13 


1.89 


4 2 11 


1.30 


6268 


637 1 


18 12 33.9 


-18 23 51 


166691 


B8 II 


8.23 


8.20 


7.92 


7.89 


7.43 





1.23 


1.61 


1.39 




43 30 


3.00 


5689 


949 1 


18 19 04.9 


-13 48 20 




B 1.5 V 


11.12 


10.69 


9.58 


9.45 


9.31 


1 


0.44 


1.05 


1.75 


2.44 


4 13 2 


3.88 


5689 


815 1 


18 19 05.6 


-13 54 50 




BIV 


9.73 


9.50 


8.71 


8.65 


8.61 





0.21 


0.90 


1.17 


4.21 


4 13 4 


2.43 


6269 


1768 1 


18 19 27.2 


-18 13 10 


168229 


B2/3 IB/nE 


8.86 


8.68 


7.90 


7.76 


7.50 





0.17 








2000 


0.98 


6265 


875 1 


18 20 22.2 


-15 48 29 


168446 


G8 V 


10.76 


9.91 


8.21 


7.86 


1.11 





0.26 








4000 


1.79 


6269 


980 1 


18 22 50.3 


-17 26 48 


313240 


B2 IV(NN) 


9.51 


9.34 


8.57 


8.44 


8.25 


1 


0.12 








2000 


1.36 


6274 


832 1 


18 27 12.0 


-18 57 13 


169805 


B2 VNE 


8.18 


8.04 


7.42 


7.30 


7.05 





0.30 








4 


1.16 


5124 


2597 1 


18 36 07.8 


-06 44 31 


171610 


K2 III 


8.67 


7.09 


4.77 


3.66 


4.42 





1.75 


0.83 


0.55 


9.55 


444 1 


0.86 


5700 


654 1 


18 39 39.8 


-11 52 42 


172252 


BO V:e 


10.12 


9.63 


7.93 


7.74 


7.49 





0.17 








3 000 


0.98 


5125 


325 1 


18 44 33.3 


-07 06 38 


173219 


BO lae 


8.04 


7.90 


7.14 


7.05 


6.80 





0.32 








4 


0.99 


1044 


195 1 


19 04 08.0 


-1-11 06 24 


230579 


B1,5:IV:NE 


9.72 


9.19 


7.74 


7.54 


7.26 


1 


0.22 








4 


1.02 


2139 


1057 1 


19 42 50.9 


-1-22 33 45 


344800 


B2V:NNE 


10.42 


10.05 


8.90 


8.78 


8.49 





0.10 


1.13 






2 10 


1.39 


2682 


3762 1 


19 59 55.2 


+31 02 34 


189687 


B3V 


4.97 


5.12 


5.41 


5.48 


5.49 





0.85 








4000 


0.83 


2679 


717 1 


20 09 58.4 


-1-35 29 45 


228041 


B0,5V:E 


9.35 


9.05 


7.87 


7.64 


7.36 





0.19 








3 000 


0.98 


2683 


1156 1 


20 13 33.0 


-1-36 19 42 


192445 


B0,5III 


7.03 


7.08 


7.13 


7.08 


6.91 





0.32 








4 


1.13 


2683 


174 1 


20 13 50.3 


-1-36 37 22 


228438 


BO(IV) 


8.77 


8.41 


7.18 


7.05 


6.76 





0.30 








4 


0.85 


2686 


693 1 


20 33 05.1 


-1-31 39 25 


195907 


B 1.5 V 


7.82 


7.79 


7.46 


7.43 


7.28 





0.19 








2 


0.96 


3171 


1117 1 


20 54 22.3 


-1-40 42 10 


199218 


B 8 V:NN 


6.59 


6.69 


6.69 


6.73 


6.70 


1 


0.26 








4000 


0.77 


3575 


553 1 


20 57 59.4 


-1-46 28 00 




M2IA 


11.66 


8.48 


2.77 


1.74 


1.29 





113.22 


119.72 


80.39 


86.49 


4444 


1.76 


3970 


673 1 


21 24 30.3 


-1-55 22 00 


204116 


BIVE 


8.13 


7.60 


6.21 


5.96 


5.63 





0.88 




0.89 




4 10 


0.88 


3194 


2107 1 


21 25 02.4 


-h44 27 06 




B1,5V:PNNE 


9.47 


8.94 


7.43 


7.21 


6.94 





0.26 








4000 


0.91 


3979 


1275 1 


21 36 59.6 


-1-58 08 24 


239712 


B2 VNNE 


9.00 


8.59 


7.54 


7.48 


7.28 





0.23 








2000 


1.12 


3986 


1670 1 


22 18 45.6 


-1-56 07 33 


211853 


BO:I:-i-WR 


9.44 


9.06 


7.78 


7.60 


7.41 





0.18 


1.18 


1.09 




3110 


1.00 


4282 


933 1 


22 56 42.6 


-1-62 37 29 


217061 


BIV 


9.51 


8.88 


7.26 


7.13 


7.03 





0.64 


1.85 


2.71 


5.70 


4244 


2.03 


4011 


773 1 


23 20 34.3 


-H58 16 39 


220116 


B0,5VPE 


9.26 


8.68 


7.59 


7.50 


7.36 





0.17 








2000 


0.90 


4285 


1070 1 


23 49 53.1 


+62 12 50 


223501 


B 1 VNE 


7.73 


7.74 


7.74 


7.77 


7.67 





0.18 








3 000 


1.28 
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APPENDIX A: OTHER MSX CATALOGUES 

The MSX satellite observed a number of regions in addition to the 
galactic plane and the IRAS gaps. These were generally regions of 
high source density that were labelled as confused by IRAS and 
included the Magellanic Clouds, several star forming regions and a 
number of extended objects such as galaxies. We have undertaken a 
search for excess stars in the IRAS gaps and all of the above regions 
except the galaxies. The results of these searches are discussed be- 
low. 

A positional correlation of the non-plane (| b | > 6°) MSX cat- 
alogue sources with the Tycho 2 star catalogue found 5898 (68%) 
MSX sources within 6" of a Tycho 2 star. The non-plane MSX cata- 
logue contains both the IRAS gap (section lAll and the Large Mag- 
ellanic Cloud (LMC) surveys. We discuss the LMC stars separately 
in section lA2l and define these to be within a 14° square (272°,- 
26°to 286°,-40°) around the LMC's galactic coordinate position. 
Of the non-plane sample 753 (12%) were within 45"of an IRAS 
source. A small fraction (20%) of these IRAS associations are LMC 
stars. 

Al IRAS gaps 

5364 MSX-Tycho 2 IRAS gap stars are identified. 66 % (3522) of 
these have spectral type information and 5215 have 2MASS coun- 
terparts. A near-mid IR colour diagram is shown in Figure |aI] Al- 
though much less reddened, the H~K K~[8] colours are similar 
to the galactic plane sample. The same excess determination pro- 
cess was applied to the IRAS gaps, and we identify 95 (87 warm 
and 8 hot) IR excess stars in this region. 22 of these excess stars 
were found to be within 45"of a IRAS source; these were found to 
be at the edges of the IRAS gap regions. 

The star with large excess is the variable star ST Pup, the other 
stars with excesses fall neatly into the groups we have previously 
discussed for the galactic plane. These stars are published in Table 
3. 

A2 Magellanic Clouds 

The LMC was found to contain 534 MSX-Tycho 2 sources, of these 
523 have 2MASS counterparts and 333 have spectral type infor- 
mation. A H—K K~\^] colour diagram of the sample is shown 
in Figure IaTI These stars are mainly G,K and M giant type stars. 
However a comparison of MSX and 2MASS observations of the 
LMC by E gan. van Dvk and Pried i2001f) indicated that the spatial 
distribution of these stars are more likely galactic foreground ob- 
jects rather than LMC stars. We identify 24 warm excess stars, 10 
of which were within 45"of a IRAS source. The LMC infrared ex- 
cess stars can be found in Table 4. 

The excess stars are nearly all supergiants, the objects with 
H~K > 0.70 are well known Be supergiants, their extreme H—K 
colours are a result of strong K band excess from the stars hot 
dust disk. Interestingly the colours of the Be supergiants overlap 
with the Herbig stars observed in the galactic plane (see Fig|4j. 
The stars with blue H—K < 0.2 colours and large excesses are 
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young massive stars. The lower excess stars are late type stars with 
circumstellar material and the well known Luminous Blue Variable 
S Dor. 

A correlation of the Small Magellanic Cloud (SMC) MSX 
mini-catalogue (containing 243 sources) with Tycho 2, yielded 55 
associations. 54 of these stars have 2MASS counterparts and 42 
have spectral types. A H—K K—^] colour diagram of the SMC 
sample is shown in Fieure lATI We identify 4 excess stars (4 warm). 
Of these none were detected by IRAS. The excess stars found in 
the SMC were known or suspected late type supergiants. The SMC 
infrared excess stars are listed in Table 5. 

A3 Star forming regions outside of the galactic plane 

Here we discuss only the regions that are outside of the galactic 
plane, and hence are not part of the galactic plane survey discussed 
earlier The follow ing star forming regio ns were imaged during the 
MSX mission (see lKraemer et alj2003l for full details): 

The Orion Nebula (A & B), the HII region S263, the IRAS 
loop G159.6-18.5 associated with the Perseus Molecular Cloud, 
the Pleiades star cluster and G300.2-16.8 an isolated cloud in 
Chamaeleon associated with IRAS 1 1538-7855. We have searched 
the MSX Point Source mini-catalogues associated with each of 
the regions for Tycho 2 and 2MASS counterparts. The searches 
found 160 MSX-Tycho2-2MASS sources in Orion, 41 in S263, 53 
in G159.6-18.5, 51 in the Pleiades and 54 in G300.2-16.8. 

We show the -R'— [8] H—K colour diagram for these star 
forming regions in Fieure lATI The majority (78 %) of the IR ex- 
cess stars in the star forming regions are found in the Orion Nebula. 
In the star forming regions we identify 51 excess stars (46 warm 4 
hot 1 cool). Of theses sources approximately 80 % were within 
45"of an IRAS point source. The excess stars tend to split into the 
Herbig stars with H—K > 0.5 and variable stars of Orion type 
H—K < 0.5. The infrared excess stars found in the star forming 
regions are in Table 6. 
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Figure Al. Near-Mid IR colour diagram of the IRAS Gaps (5215 stars) (left),the LMC □ (523 sources) SMC * (54 sources) (centre) and the non-galactic 
plane star forming regions observed by MSX (right). In the right hand panel, the different regions are indicated by the following symbols: x G159.6-18.5, A 
G300.2-16.8, □ Orion Nebula (A and B), o Pleiades and a * for S263. The excess cutoff is shown in all the diagrams by a broken line. 



